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Abstract 

The quantum master equation is introduced for the density matrix representing Bogoliubov-BCS 
quasiparticles. A constraint to relate the loss and gain factors is taken into account to preserve the 
form of the density matrix. Such an equation can be reduced to the semiclassical equation, and 
can be extended for the coexistence of different order parameters. 

PACS numbers: 
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I. INTRODUCTION 



16| 



Considerable efforts have been made to develop quantum kinetic approaches, [l 
Although semiclassical kinetic approaches can be used to model many irreversible processes, 
quantum corrections should be taken into account when we consider the systems of the nano- 
scale. Different master equations have been introduced for such corrections by including 
quantum relaxation terms in addition to the Liouville term. [l-9| The Markoff master 
equation is successful in quantum optics, 2| and the equation of Lindblad form js, 4| is 
derived by considering suitable assumptions. It is also discussed in the literature how to 
construct the quantum master equations to include Fermi and Bose properties, js, 6| The 
nonlinear relaxation terms have been introduced in Refs. [sl and {t] for non- interacting 
identical fermions, and we can see the equivalence after some calculations. In addition to 
the formal derivation, we can obtain such nonlinear terms intuitively by considering both 
the lifetimes of particles and holes based on the conservation of the number of particles in 
each transition, js, 8| Here holes are vacancies of any orbitals. The lifetimes of particles and 
holes describe the loss and gain of particles, respectively, and the equations introduced in 



7| for fermions are symmetric with respect to particles and holes. 



Refs. 

Because the nonlinear relaxation terms in Refs. |6| and fl are for non-interacting 
fermions, it is natural that they can become invalid when many-body effects are impor- 
tant. [9] Such effects have been taken into account under significant two-body interactions, 
which may induce the formation of Bogoliubov-BCS quasiparticles 17H23l|. Such fermionic 
quasiparticles have been successfully introduced to understand superconductors. While the 
number of particles is conserved in the conventional noninteracting models, the conservation 
is not necessary for quasiparticles. Semiclassical master equations have been proposed to 



model Bogoliubov-BCS quasiparticles when the wave properties can be ignored. 2J, |25| On 



the other hand, WKB approximation 



24 



26| based on Bogoliubov-de Gennes equation 



27| 



provides an analytic way to understand the wave behaviors. [28] To describe the quantum 
wave properties when the relaxation effects are incorporated, it is of the fundamental interest 
to extend the quantum master equations to Bogoliubov-BCS quasiparticles. 

In this paper, the quantum master equation is extended for the fermionic density matrix 
representing Bogoliubov-BCS quasiparticles. Without the loss of generality, we focus on the 
quasiparticles due to the pairings of superconducting electrons in the coordinate space. In 
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section II, we discuss how to construct the loss and gain factors for such a density matrix 



23| 



based on the general form introduced in Ref. [g]. Different types of order parameters 
can be taken into account in the extended equation, as shown in section III. The conclusion 
is made in section IV. 



II. A QUANTUM MASTER EQUATION FOR BOGOLIUBOV-BCS QUASIPAR- 
TICLES 

The master equation of the following form 

= i[m.m] - {mm) + {/ - m. ^'m (i) 

has been introduced in Ref. for identical fermions. Here f){t) and H{t) represent the 
density matrix and Hamiltonian for the considered fermions, / is the identity operator, 
and we denote [Oi,02] and {Oi,02} as the commutator and anticommutator for any two 
operators Oi and O2. In this paper, we take the reduced Planck constant h = 1. We require 
that f (t) and f '(t) are both positive self-adjoint operators such that the above equation can 
preserve the positivity and Pauli's exclusion principle under suitable assumptions. |i8|] The 
operator / — p{t) represents the density matrix for the corresponding holes, and the last two 
terms in Eq. (1) describe the lifetimes of particles and holes for the^ irreversible effects due 

'or bosons, we just 



7| that the form 



to the relaxation and/or excitation. [6|] To extend Eq. (1) as that js, 6] 
need to replace the last term by {/ + p(t), V'{t)}. We can see from Refs. 
of Eq. (1) can be derived by different ways. 

To see the meanings of the last two terms in Eq. (1) explicitly, consider the noninteracting 
spin-unresolved electrons in a finite cube of the volume V under the periodic boundary 
condition. In addition, assume that the Hamiltonian can be approximated as the time- 
independent operator Hg = Xl^kCko-Cko- around a specific time ti. Here Ckf (cki) denotes the 
annihilator for the spin-up (spin-down) electron with the plane wave tt(r; k) = p^exp(2k-r) 
as the wavefunction, k is quantized because of the boundary condition, and the real number 
represents the eigenenergy of the spin-unresolved orbital corresponding to Cko-. Because 
of the unresolved spin-splitting, we can set (cj^^Ckv) = {c]^-a^k'-cT')Sa,a' and do not need to 
consider the spin-orientation. Here we take (A) as the expectation value of any (many-body) 
operator A with respect to the total density matrix, which describes both the reservoir and 
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the considered system. When the transitions are between the eigenorbitals of Hg around a 
specific time ti, we shall set the operators t(t) and t'(t) in Eq. (1) as 

f,(t) = 5^a,(k;tK(r;k)<(r';k) and t'^{t) = /3,(k; tK(r; k)<(r'; k) (2) 

k k 

in the coordinate space when t ~ ti. Here as(k;t) and /3s(k;t) are the positive coefficients 
for the loss and gain in orbital k as t ~ ti, and for any complex number z we denote z* 
as its complex conjugate in this paper. Let Fs(k;t) = {c]^^Cka) as the number of electrons 
annihilated by Cko- at time t. We have from Eqs. (1) and (2) 
d 



dt 



F,(k; t) = -2a,(k; t)F,(k; t) + 2/3,(k; - F,(k; t)). 



(3) 



The function Fs{k]t) follows < Fs(k;t) < 1 because electrons are fermions, and the first 
and second terms at the right hand side of the above equation is to decrease and increase 
the occupation number in orbital k. In Eq. (1), therefore, the operators f (t) and f'(t) 
responsible for these two terms can result in the loss and gain of electrons, respectively. The 
above equation is just the nonlinear semiclassical master equation 29|] if 

a,(k; t) = ^ ^(1 - F,(k'; t)) and /3.(k; t) = J2 "^m'; t)- (4) 

k' k' 

Here Wkk' denotes the positive coefficient for the transition from k' to k if k 7^ k', and we 
take Wkk = for all k. Because the semiclssical master equation is not linear, it is natural 
that the last two terms in Eq. (1) is nonlinear, j^, Q] 

To extend Eq. (1) for Bogoliubov-BCS quasiparticles in superconudctors, we note that 

/i 



the quasiparticle wavefunction ff, 



is the direct sum of an electron and a hole 



/2 



components, where /i and /2 are defined in the coordinate space. For convenience, we 
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27| 



(5) 



denote = {f^ , f2)- It is known that the effective Hamiltonian is of the form 

where Hs{t) denotes the single-electron Hamiltonian, /i is the chemical potential, and Ks{t) 
represents the pairing field. The two signs in ± have been introduced for different BCS 
models. The effective density matrix is 211] 

Psit) Asit) 

±A:(t) is-m 



pb{t) 



(6) 
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where Ps{t) denotes the one-body density matrix for electrons, As{t) is the pairing tensor, 
and is = ig is the identity operator for single-electron orbital. The operators Hb{t) and 
Pb{t) follow 



SbPi{t)si=h-Pb{t) 



(7) 



-Hh{t). 



(8) 



Here h 







Is 
L 



is the identity operator for Bogoliubov-BCS quasiparticles, and Sb 



We can see from the direct calculations that the Hamiltonian and density 

T/s " 

matrix satisfying the above two equations can be wriiten as the forms given by Eqs. (5) and 
(6). To model Bogoliubov-BCS quasiparticles by Eq. (1), we shall take pb{t) and Hb{t) as 
the density matrix and Hamiltonian in such an euqation. Let tb{t) and f^t) as f (t) and 
r'(t) for the lifetimes of quasiparticles and quasiholes. It will be shown in this section that 
we shall introduce the constraint 



Sbti{t)si = n(t) 



(9) 



to preserve Eq. (6). The lifetimes of Bogoliubov-BCS quasiparticles and quasiholes are 
related under the above equation, which is reasonable because we can obtain the density 
matrix for quasiholes from that for quasiparticles by Eq. (7). 

To see why we need to introduce the above constraint, first we consider a s-wave pairing 



case where the Bogoliubov transformation 



2Q|, 



25| 



7kt = ^^kCkt + 4clk; 



= -44t + ^kC_k; 



(10) 



yields the annihilators for the quasiparticle with the excited energy ^k > around a specific 
time ti- Here the complex numbers zy_ and z'-^^ denote the transformation coefficients. Let 
Ub(r; k+) and U5(r; k— ) = S';,u^(r;k+) be the positive- and negative-energy orbitals corre- 
sponding to the quasiparticles annihilated by 7^0-, and assume that the Hamiltonian can be 
approximated as the unperturbed one Hb = Xlk '^k[uf,(r; k+)u^(r'; k+)— Ub(r; k— )u||(r'; k— )] 
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as t ~ ti. The number of quasiparticles in orbital k is Fb(k] t) = tr[/3fe(t)ufe(r; k+)u^(r'; k+)]. 
We also have Fb(k; t) = I — tr[pb{t)ub{r; k— )u^(r; k— )] because 

tr[pb{t)ub{r; k+)u;(r'; k+)] = 1 - tr[pb{t)nb{r; k-)uj(r; k-)] (11) 

from Eq. (7). Let 2Q;f,(k, t) and 2/35(k, t) be the loss and gain rate of the quasiparticles 
in orbital k when t ~ ti. The loss (gain) in orbital k not only decreases (increases) the 
occupation number in Ub{r; k+), but also increases (decreases) of the occupation number in 
Uft(r; k— ) because of the above equation. By checking the loss and gain rates, we shall set 

hit) = J2 «b(k; t)ubir; k+)u;(r'; k+) + f^)(t) (12) 

k 

hit) = 5^/3,(k;t)u,(r;k+)ul(r';k+) + f[^-)(f) (13) 

k 

with the two operators 



^ b 



Sb Mk- t)ub(r; k+)u;(r'; k+) j Sl = /3,(k, t)u,(r; k-)u;(r'; k-) (14) 

h''^ = Sb Mk, t)u,(r; k+)ui(r'; k+) j §1 = ^ a,(k, t)u,(r; k-)ui(r'; k-). (15) 

In comparison with Eq. (2), the first terms at the right hands of Eqs. (12) and (13) can 
be taken as the extensions of f ^ and f to BCS models by Ms(r; k) — )• Ufe(r; k+), as(k; t) — )• 
af,(k;t), and /3s(k;t) — )■ /3bik;t). On the other hand, f^'"-* and f[^~-* are to relate the 
occupation numbers in negative- and positive-energy orbitals. We can see that Eqs. (12) 

q. (9) by direct calculations. 
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should be incorporated in addition to 



and (13) yield operators Tbit) and r'{,(t) following 
It is known that pairing creations/annihilations 
quasiparticle transitions. By checking the loss and gain rates, we can obtain the coefficients 

(k;t) = ^ E4'.!(l - ^^(k'; t)) + ^ $^42F,(k'; t) (16) 



k' k' 



/3,(k;t) = ^$^42^^(k';t) + -^.(k';t)) (17) 



2 ^ ' ' 2 

k' k 



(T) (A) (C) 

for the relaxation term due to the electron-phonon interaction. Here wj^,^,, uj\^y/,, and u^^J 
represent positive real parameters for the quasiparticle transition from k' to k, pairing 



annihilation for quasiparticles in k and k', and pairing creation for quasiparticles in k and 
k', respectively. Under Eqs. (12)-(17), the quantum master equation can be reduced as 

^F6(k, t) = -2a,{k; t)F,{k, t) + 2/3,(k; - F,{k, t)) (18) 

k' k' 

+ J2^T'Fb{k';m - F,{k,t)) + - F,{k';t)){l - F,{k,t)). 

k' k' 

The above equation, in fact, is just the semiclassical master equation for Bogoliubov-BCS 
quasiparticles when the relaxation term is due to the electron-phonon interaction. 24 1 
In general, we can include the spin orientation and extend Eq. (11) as 

tr(p,(t)f4) = 1 - tr[Mt){s,f:){s,f:y] (19) 

for any normalized quasiparticle wavefunction f;,. If f ^(t) (f ^t)) induces the decrease (in- 
crease) of the number of quasiparticles in f;,, t[(t) {tb{t)) should induce the increase (de- 
crease) of the occupation number in Sbi^ based on the above equation. Therefore, we shall 
use Eq. (9) to relate Tb{t) and rj,(t). When Eq. (8) is valid, the form given by Eq. (6) 
is preserved under Eq. (9) if the time evolution of the density matrix follows Eq. (1). To 
see this, assume that a fermionic density matrix pb{t) follows Eq. (1) with H(t) = Hb{t), 
i = ib, f (t) = tb{t), and t'(t) = t'bit)- (Here an operator O is fermionic iff the inner 
product {a,Oa) for any normalized ket a is a real number between and 1.) In addition, 
assume that SbPl(ti)Sl = h — pb{U) such that Eq. (7) holds true at the initial time ti, and 
let f)'^{t) = Sbih- P*b(t))Sl = ib- SbPl{t)Sl when t > U. It is easy to see that p^tj) = pbiU), 
and p'b{t) is also fermionic. The time derivative of p'b{t) follows 

|pUt) = iSb[pl{t),Hl{t)]Sl + Sb{pl{t)Jl{t)}Sl (20) 

-Sb{ib-mx:msi 
= Hb{t)] + [ib - m. f ut)} - {mMm 

from Eqs. (8) and (9). Hence pb{t) is also a fermionic density matrix following Eq. (1) with 
the initial condition the same as that for Pb{t)- Because the uniqueness of the solution to 
Eq. (1) is expected under suitable assumptions, [s] we have pb{t) = p[{t) = Sb{h — Pbi'^))Sl 
for any time t > to- The equality implies Eq. (7), under which pb{t) is of the form given by 
Eq. (6). Therefore, the form of Pb{t) is preserved under Eq. (1) when Eqs. (8) and (9) are 
valid. 
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III. THE MASTER EQUATION FOR THE COEXISTENCE OF DIFFERENT OR- 
DER PARAMETERS 



In addition to superconducting electrons, other identical fermions may also form 
Bogoliubov-BCS quasiparticles by pairing. Equation (6) can be used for different Fermi 

systems if we replace ps(t) and A(t) by the density matrix and pairing tensor for the corre- 
sponding fermions. For fermionic Bogoliubov-BCS quasiparticles, we can substitute pb{t) for 
Ps{t) at the right hand side of Eq. (6) to extend such type of quasiparticles by introducing 
the density matrix 



Pb{t) Ah{t) 

±Ai{t) h-m 



(21) 



Here Abit) is the pairing tensor for the quasiparticles described by pb{t), and pBit) represents 
the extended Bogoliubov-BCS quasiparticles. Because pbit) and lb — plit) are 2x2 matrices, 
A(,(i) is also a 2 x 2 matrix and pB{t) is a 4 x 4 matrix. It will be shown in this section that 
different order parameters can be incorporated by introducing psit), and the following two 
constraints 



SBPms^B-iB-Mt) 



(22) 



ABPB{t)^B = PB{t) 



for pB{t) will be introduced based on Eq. (7). Here the 4x4 matrices Sb 



h 
h 



, and Ab — 



Sb 
^5-6 



(23) 



Sb 
Sb 



. The operator Ib is just the identity matrix for 



the particles described by ps(t). To model pB{t) by Eq. (1), let HB{t), fs(t), and fs(t) as 
the corresponding H{t), r(t), and t"{t) in such an equation. It will be also shown in this 
section that we shall introduce 



5s^;(t)4 = -HB{t) 



(24) 



SBtut)sl, = t'^it) 



(25) 
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for the validity of Eq. (22), and introduce 

ABHB{t)A^B = HB{t), 



(26) 



ABrBit)Al^rB{t) 



(27) 



for Eq. (23). 

Because Eq. (21) is obtained from Eq. (6) by taking pi,{t) as Ps{t), we can see that the 
constraint 



s'BP*B{t)s'^ ^ Ib - Mt) 



(28) 



with the 4x4 matrix 5'r 



is equivalent to Eq. (21) from the equivalence 



h 
Th 

between Eqs. (6) and (7). In addition, the following two constraints 



S'BH*B{t)Si - -HB{t) 



(29) 



(30) 



should be introduced for Eq. (28) just as Eqs. (8) and (9) are introduced for Eq. (7). To 

use Pb{t) to represent the quasiparticlcs resulting from electron pairing, however, we shall 
introduce additional constraint for Pb{t) to follow Eq. (6). 

We can see from the last section that Eq. (6) is preserved iff Eq. (7) holds true, and Eqs. 

Sb 



(8) and (9) are important to the vahdity of Eq. (7). The 4x4 matrix Sb 



Sb 



is the 



natural correspondence to the 2x2 matrix Sb, so it is reasonable to extend Eqs. (7)-(9) as 
Eqs. (22), (24), and (25). Prom Eq. (22), Ab{t) is of the form 



Kit) m 



(31) 



The matrix Ab equals the product SbS'b- Therefore, 

ABPBmi = Sb{S'bpW)S'^)SI = hit) 



(32) 



and we can obtain Eq. (23) from Eqs. (22) and (28). In fact, Eq. (23) is equivalent to Eq. 
(28) under Eq. (22), so we only need to consider Eqs. (22) and (23) for the form of hit) 
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and Pb{t)- By checking ^^^^(t)^^ and ABTB{'t)A^B^ can see that Eqs. (26) and (27) 
are equivalent to Eqs. (29) and (30), respectively, if Eqs. (24) and (25) hold true. So we 
just need to take Eqs. (24)- (27) as the constraints on the quantum master equation. By 
considering the time-derivative of Sb{Ib — P*Bi^))^B ^si^B — P*Bi^))^w '^^^ prove 
that the form of pB{t) is preserved under these constraints. 
From Eqs. (6), (21), and (31), we can rewrite Pb(^) by 



hit) 



(33) 



(34) 



Ps{t) A,(t) A',(t) 6s{t) 

±A:(t) is-m ±6:it) -A-(t) 
±A':it) ±5:it) is-m -kit) 

Assume that be a normalized eigenfunction of psit) at a specific time t = ti. Because 
PB(ti) is a 4 X 4 matrix, includes four component and we can write 

fhl 
fh2 
\fe2j 

By checking the contribution of f^f^ to pB(ti), we can see that the components feifei 
/e2/e*2 ^irc incorporated in the first and fourth diagonal terms, both of which are just p.s{ti). 
On the other hand, fhifhi fh2fh2 incorporated in the second and third diagonal 
terms, both of which are Ig — p*s{t)- Because Ps{ti) and Ig — Ps{ti) represent the density 
matrices for electrons and holes at ti, we shall take fei and /e2 as electron components and 
take fhl and fh2 as hole components. Therefore, contain two electron components and two 
hole components just as the four-component wavefunctions introduced for the coexistence 

I. The coexistence may 



of the superconducting and antiferromagnetic orders in Ref. 
reveal the key to understand the high-temperature superconductors. The operator 6s(t) in 
Eq. (33) can correspond to the antiferromagnetic order while A'^(t) and As{t) can be taken 
as the superconducting orders, and the extended master equation for pB^t) could be used 
to model the nonequilibrium phenomena when there are different order parameters. In the 
stationary case, it may be convenient to remove the degeneracy indicated by in Eq. (23). 

For the BCS-type quasiparticles described by PB{t), thejdensity matrix for the corre- 
sponding quasiholes is Ib — PB{t)- Just as mentioned in Ref. 



, we can substitute Pb(0 for 
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Ps{t) at the right hand side of Eq. (6) to the further extension. A chain of density matrices 
can be constructed based on Eq. (6). [s] On the other hand, we can extend the Bogohubov- 
BCS density matrix by considering the quasiparticles with multiple electron components 
and the coupling between them and the corresponding quasiholes. [s^ The orbial in Eq. 
(34), in fact, can be interpreted as coupling between the antiferromagnetic-like quasiparticles 

"^^^ I and the corresponding quasiholes j '^^^ j . There is no upper limit to the number 

fe2 J y A2 y 

of components in principal, and it may be convenient to introduce the extra dimensions for 
the quasiparticles with the infinite components. 



IV. CONCLUSION 

A quantum master equation is obtained for the density matrix representing Bogoliubov- 
BCS quasiparticles. Such an equation can be reduced to the semiclassical equation, and can 
be extended for the coexistence of different order parameters. 
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